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Abstract 

We explore a Z' boson with family-nonuniversal couplings to charged leptons. The general effect of 
Z-Z' mixing, of both kinetic and mass types, is included in the analysis. Adopting a model-independent 
approach, we perform a comprehensive study of constraints on the leptonic Z' couplings from currently 
available experimental data on a number of flavor-conserving and flavor-changing transitions. Detailed 
comparisons are made to extract the most stringent bounds on the leptonic couplings. Such information 
is fed into predictions of various processes that may be experimentally probed in the near future. 
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I. INTRODUCTION 



Recent anomalous measurements of a number of observables at the Fermilab Tevatron, such 
as the forward-backward asymmetry in top-quark pair production [1], the like-sign dimuon charge 
asymmetry in semileptonic 6-hadron decays [2], and the invariant mass distribution of jet pairs 
produced in association with a W boson [3], give us possible hints on physics beyond the stan- 
dard model (SM). One of the candidates that have been proposed to explain these anomalies is 
a massive spin-one electrically neutral gauge particle, the Z' boson, which may be associated with 
an additional Abelian gauge symmetry, U(l)', that is broken at around the TeV scale and has 
a mass of —150 GeV [4-6]. Moreover, the desired Z' boson would need to have sufficiently sizable 
flavor-changing neutral-current (FCNC) interactions in the quark sector. 

One way to induce Z'-mediated FCNC's is to introduce exotic fermions having U(l)' charges 
different from those of the SM fermions [8], as occurs in models with the E 6 grand unified group. 
In this case, mixing of the right-handed ordinary and exotic quarks, all SXJ(2) L singlets, gives rise 
to FCNC's mediated by a heavy Z' or due to small Z-Z' mixing. Another possibility involves 
family-nonuniversal interactions of the Z'. In string-inspired model building, it is natural for at 
least one of the gauge bosons of the extra U(l) groups to possess family-nonuniversal couplings 
to ordinary fermions [9]. In this scenario, the FCNC couplings appear when one transforms the 
SM fermions into their mass eigenstates, without the necessity to introduce new fermion states. 
Furthermore, both left- and right-handed fermions can have significant flavor-violating interactions 
with the Z', as well as small family-nondiagonal couplings to the Z boson caused by Z-Z' mixing. 

In fact, Z 1 models with tree-level quark FCNC's have been studied extensively in low-energy 
flavor physics phenomena, such as neutral meson (K, D, or B) mixing, B meson decays involving 
the b — > s transition in particular, and single top production [10-12]. In principle, one can consider 
the possibility of FCNC's in the lepton sector as well [10]. In this work, we focus on family- 
nonuniversal interactions of the Z' with the charged leptons and explore constraints on its relevant 
couplings from various experiments involving only leptons in the initial and final states. Such 
processes suffer less from QCD corrections and hadronic uncertainties than the above-mentioned 
hadronic systems. We assume that the Z' boson arises from an extra U(l) gauge symmetry, but 
otherwise adopt a model-independent approach. We take into account the effect of Z-Z' mixing, of 
both kinetic and mass types, which modifies theoretical predictions of the electroweak p parameter 
and various Z-pole observables. Due to the family nonuniversality, such a Z' boson would feature 
flavor-changing leptonic couplings, as would also the Z boson through the mixing. We therefore 
examine a number of flavor- conserving and flavor-changing processes to evaluate constraints on the 
leptonic Z' couplings. 

This paper is organized as follows. We present the interactions of the Z' boson with the charged 
leptons in Section II. The p parameter from global electroweak fits is used to determine the allowed 
mixing angle between the Z and Z' . In Section III, we study constraints on the flavor- conserving 
couplings of the Z' . The pertinent observables include those in leptonic Z decays from the Z-pole 
data and the cross sections of e + e~ collisions into lepton-antilepton pairs measured at LEPII. 
We separate the analysis of the flavor-changing couplings into two parts. The constraints from 
transitions generated by tree-level diagrams are treated in Section IV. We place upper bounds on 
the couplings from the rates of flavor-violating Z — >■ W decays, p — > 3e, several flavor-violating 
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t decays into 3 leptons, muonium-antimuonium conversion, as well as the cross sections of flavor- 
changing annihilations e + e~ — > 11'. The constraints from processes induced by loop diagrams 
are given in Section V. The considered processes or observables are the flavor- changing radiative 
lepton decays / — > the anomalous magnetic moments of leptons, and their electric dipole 
moments. We will make use of the existing experimental information on all these transitions, 
including new measurements from the BaBar, Belle, and MEG Collaborations [13-15]. Based on the 
allowed coupling ranges, we make predictions for various flavor-conserving and -violating processes 
in Section VI. These predictions can serve to help guide experimentalists in future searches for Z' 
signals. Our findings are summarized in Section VII. 



II. INTERACTIONS 

The mass Lagrangian for the interaction eigenstates Z and Z' of the massive neutral gauge bosons 
after electroweak symmetry breaking, which leaves the photon massless, can be expressed as 




(1) 

where M zz , denote the masses of the gauge bosons and A represents the mixing between them. 
As discussed in Appendix A, which has some more details on the notation we adopt, A contains 
both possible kinetic- and mass-mixing contributions, and in the presence of kinetic mixing the 
parameter M z , is not identical to the orignal mass of the U(l)' gauge boson [see Eq. (A7)]. 

The squared-mass matrix in £ m can be diagonalized using [16] 

Z\ ( cos£ -sin A ( Z\ 2A 

with its eigenvalues being 

m%,z' = \{M 2 Z + Ml) =F i V / (Ml-Mf,) 2 + 4A2 . (3) 

One can then derive 

(m|,-M|)tan 2 e = M\-m\. (4) 

The Lagrangian describing the interactions of Z and Z' with the charged leptons is 

Ant = ~9zJz^\ — 9z>Jz>Z'\ , (5) 

and the currents are given by 

9z Jl = (g L P L + 9rPr) i , 9z> Jz> = {sf L P L + sTrPr) L (6) 

where I — (e fi f) T contains the interaction eigenstates of the leptons, P LR = |(1 =F 75), and 
the coupling constants g L R are family universal, whereas the Z' couplings are not assumed to be 
family universal according to 

g' L = diag(L' e ,L^,L;), g' R = diag(^, R'^ R' T ) , (7) 
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with the parameters L' and R' e fl T being generally different from one another. The Hermiticity 
of C int requires these coupling constants to be real. The interaction eigenstates in £ are related to 
the mass eigenstates in £ = (e \i r) T by 1 

i L = P L £ = V L £ L , 4 = Pj = V R £ R , (8) 

where V L R are unitary matrices which diagonalize the lepton mass matrix Mi in the Yukawa 
Lagrangian, diag(m e , m , m T ) = V'[M ( V R . 

In terms of the mass eigenstates, Z, Z' , and £, we can then write 

Ant = -h X [{g L + B L sin£)P L + (g R cos£ + ^ sin $P R ]£Z X 

- h X [(-9 L sine + cose)P L + sin£ + 5 R cos^)P R ]£Z' x 

= -IffffiPL + Pr'Pr)^ Zx ~ W^Pl + ^Pr)^ Z' x , (9) 

where B L = V[ g' L V L and B R = V R g' R V R are generally nondiagonal 3x3 matrices, summation over 
i,j = 1,2,3 is implied, £1,2,3 = e, /x, r, and 

= = ^-c c (/ c + Sc (S c ) ., 6^ = -^. ac (/ c + Cc (S c ) . (10) 

for C = L or i?, with = cos£ and = sin£. One can see from Eq. (9) that the presence of 
nonzero off-diagonal elements of B LR , due to the nonuniversality of the diagonal elements of g' L R 
and to the charged-lepton mixing, gives rise to flavor-changing couplings of the Z' to the leptons 
at tree level. Furthermore, Z-Z' mixing introduces not only family nonuniversality, but also flavor 
violation into the tree-level interactions of the Z. 



Now, it follows from Eq. (10) that 



= + , (11) 



where = tan£. Therefore the couplings of Z and Z' to £ i £j are directly related once the mixing 
angle £ is specified. Employing the electroweak data, one can fix £ if the Z' mass is given. We 
achieve this by means of the p parameter, which in the Particle Data Group (PDG) convention [17] 
encodes the effects of new physics if it deviates from the SM expectation p$ M = m^/ (c^Mf ) = 1, 
where c w is the cosine of the Weinberg angle 6w Since Z-Z' mixing alters the Z mass, as indicated 
in Eq. (4), and hence causes p to shift from unity, we have 



Po 



4 m z c l M z 



ml, - Ml 2 .. 

1 ~tp3 — tan 2 £ 

Mf 



+ (12) 

Illy 



The value p = 1. 00081^00^7 \\^\ resulting from the PDG global electroweak fit then translates for 

0.008 < \i\ < 0.038 . (13) 



m z , = 150 GeV into 



Throughout the paper we make a distinction between t and I, with the former referring to the triplet of charged 
leptons and the latter to individual charged leptons in general. 
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More generally, Fig. 1 shows the corresponding limits of |tan£| for 100 GeV < m z , < 2 TeV, which 
is the range of interest in this paper. It is then straightforward to realize that for this mass range 



z tan 2 £ < 1 . 



(14) 



The plot also illustrates that |tan£| oc l/m z , as m z , becomes large, which reflects the relation 

(15) 



tan C| 



m z p - 1 



m 7J 



valid for m 2 z , 3> m 2 z and derived from Eq. (12). 



55 0.010 
8 0.005 
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FIG. 1: Values of |tan£| for 100 GeV < m z , < 2 TeV corresponding to the p range from the electroweak 
global fit. 



III. FLAVOR-CONSERVING COUPLINGS OF Z' 



With £ known, one can evaluate b l [ R from the Z-pole data. The amplitude for the Z decay into 
a charged-lepton pair is 



M 



(16) 



in the parametrization of Eq. (9). This leads to the forward-backward asymmetry at the Z pole 
and decay rate 



.4 



(o,o 

FB 



Z-tl+l- 



where 



\M 



z->i+i- 



\Jm 2 z — imf 
167r m 2 z 



W) 2 +(/3») 2 (m|-mf) +4m^M- (18) 



M 



z->i+i- 1 ' 



(17) 
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These formulas along with Eq. (11) allow us to extract b l [ R for each value of £ from the A l and 
r z _ >/+i _ measurements [18], 

Af p = 0.1515 ± 0.0019 , A e *v = 0.142 ±0.015, Af p = 0.143 ±0.004, 

r z^ e +e- = 83.91 ±0.12 MeV , r^ +M _ = 83.99 ± 0.18 MeV , (19) 
rt x P + _ = 84.08 ± 0.22 MeV , 

after # L H are fixed to their SM predictions [17] 

A su = A su = a sm = 0.1475 ± 0.0010 , 
r|^ e+e _ = r|^ M+/i _ = 84.00 ± 0.06 MeV , r|^ T+T _ = 83.82 ± 0.06 MeV . (20) 

We can reproduce all these SM numbers within their errors using Eqs. (17) and (18) with and 
/3 l R replaced, respectively, by the effective couplings 

g f = -0.1996 , gf = 0.1721 . (21) 

For comparison, their tree-level values are g L = g(s^ — l/2J/c w ~ —0.2002 and g R = gs^/c w ~ 
0.1722 if = 0.23116 [18]. We will ignore the uncertainties in gf^ R compared to the greater 
relative uncertainties in the data. 

Applying f3 l [ R = g LtR /c^ + t^b l [ R from Eq. (11) in the A l and F z _^ l+l _ formulas above, with 
g L R = g c ^ R from Eq. (21) and a specific value of £, one can then obtain the allowed ranges of b l [ R 
Thus in the m z , = 150 GeV case, for which 0.008 < |£| < 0.038, the results for £ > are 

-0.071 < bf < 0.006 , -0.11 < &|f < -0.009 , 

-0.13 < b^ < 0.25 , -0.15 < V£ < 0.27 , 
-0.070 < b T L T < 0.083 , -0.002 < b T R T < 0.16 . (22) 

Flipping the sign of £ would also flip the signs of these b u L R numbers, and the same statement 
applies to the rest of our analysis. The plots in Fig. 2 illustrate the allowed b\ R regions for the 
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FIG. 2: Values of b l [ R for m z , = 150 GeV and mixing angle £ = 0.008 (lighter colors), 0.038 (darker 
colors), as described in the text, subject to constraints from A x and T z _^i+i- data. 
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lower and upper limits of the £ range in this case, £ = 0.008 (lighter colors) and £ = 0.038 (darker 
colors). The green regions satisfy the A l constraints, red the Tz^i+i- constraints, and blue both 
of them. The upper and lower limits of the b\ R ranges in Eq. (22) are visible on the plots. 

Since Z'- mediated diagrams can also affect the collision e + e _ — > , it is important to consider 
the relevant data to see if they offer additional restraints on the Z' couplings. Here we will employ 
LEP-II measurements at various center-of-mass energies above the Z pole, from 130 to 207 GeV [19]. 
The amplitude for this process if I ^ e is 

elhne lv e If (^P L + ^P R ) I e lv (p?P L + p%P R ) e 

^e-^R ~ + m 2 z _ s 

h" {b l [P L + b ll R P R ) I ery v (bfP L + b%P R ) e 

where e p > is the proton's electric charge, s = {p e + +p e -) 2 , and we have assumed that s 
is not near m 2 zz ,. There are also contributions to this amplitude from t-channel diagrams with 
flavor-changing couplings (5f R and bf R) but we will neglect their effects in order to explore the 
largest impact of the Z' flavor-conserving couplings under the assumption that there is no un- 
natural cancellation between the two sets of contributions. Moreover, as we demonstrate below, 
the magnitudes of the latter couplings have looser upper-limits than their flavor- changing coun- 
terparts by at least a few times. Complete expressions for the cross-section <r(e + e~ — > l + l~) and 
forward-backward asymmetry Afb, including finite-width effects, are collected in Appendix B. Nu- 
merically, we adopt the couplings in Eq. (21) and the effective value a = 1/132, which in the 
absence of the Z' lead to a and Afb numbers differing by no more than 2 percent from the cor- 
responding SM predictions quoted in the LEP-II report [19]. Since we consider m z , = 150 GeV 
and larger masses from 0.5 to 2 TeV, in determining the g l [ R bounds we take the LEP-II data 
belonging to y/s = 136, 161, 205, 207 GeV for definiteness. 

We find that incorporating the LEP-II information brings about significant modifications to 
some of the results in Eq. (22). The allowed values of the couplings for m z , = 150 GeV now 
become 

-0.071 < bf < 0.006 , -0.10 < b e l < -0.009 , 
-0.033 < bf < 0.080 , -0.029 < 6^ < 0.095 , 
-0.070 < b T L T < 0.024 , < b T R T < 0.083 . (24) 

We have also explored the situations for higher masses up to m z , = 2 TeV. The inclusion of the 
LEP-II data again provide important extra restrictions on the couplings. For the representative 
values m Z i = 0.5 - 2 TeV, the allowed ranges associated with each flavor turn out to be roughly 
proportional to the m z , values, namely 
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where the numbers are in units of GeV -1 . It is worth noting that the proportionality of these 
ranges to the Z' mass for m z , ^> m z is a reflection of the |tan£| oc l/m z , behavior in Eq. (15) 
which starts to manifest itself when m z , exceeds 200 GeV or so, as can be seen in Fig. 1. We also 
note that for m z , > 2 TeV the limits in Eq. (25) accommodate couplings which may exceed order 
one in magnitude and hence the perturbativity limit. Nevertheless, as the errors in p decrease 
with increasingly better precision in future data, the bounds on b l [ R will likely become stronger. 

Before proceeding to the flavor-changing sector, a few comments regarding the case of no Z-Z' 
mixing, £ = 0, are in order. If one goes beyond the one-sigma range of the p -parameter from 
the global electroweak fit, so that the lower bound of p reaches zero, then the lower bound of |£| 
will also reach zero. In that limit R — > g L R , and therefore the Z-pole data on A l and V Z -^i+i- 
no longer offer restrictions on b\ R through the tree-level relations in Eqs. (17) and (18). At the 
one-loop level, however, Z'-mediated radiative corrections contribute to the Zl + l~ vertex, and so 
these observables can still constrain the couplings [20]. With the formulas given in Ref. [20] for the 
Z' loop contribution, we find that the upper limits on the coupling-to-mass ratios, b/m z ,, are of 
order 1 to 2 per mill for our m z , range of interest and thus higher than their counterparts in the 
presence of mixing. Without the mixing, Z'-mediated diagrams can still affect e + e~ — > l + l~ at tree 
level, as Eq. (23) indicates. The expressions for the cross section and forward-backward asymmetry 
in Appendix B suggest, however, that the LEP-II data would not impose additional restrictions in 
this case. 



IV. CONSTRAINTS FROM TREE-LEVEL FLAVOR-CHANGING PROCESSES 

A. Z -»■ e ± At = F, Z -»■ e ± r^, and Z -> 

As C int in Eq. (9) shows, the Z can have tree-level flavor-violating interactions with leptons in 
the presence of Z-Z' mixing. Accordingly, the amplitude of the decay Z — > if for V ^ I is 

M z ^f = hx(P l lP L + PiP R )l'ez, (26) 
where (3 l [ R = b l [ R from Eq. (11). The rate of this transition is then 

+ Am l m l/ Re(bfb%)Y (27) 

where p l is the three-momentum of I in the Z rest frame. These decays, like all other lepton- 
flavor-violating ones, have not yet been observed. But there is some experimental information 
available on the branching ratios: B(Z — > e ± /i =F ) < 1.7 x 10 -6 , B(Z — > e ± r =F ) < 9.8 x 10~ 6 , and 
B[Z — \ /U ± r =F ) < 1.2 x 10~ 5 [18], each of the numbers being the sum of contributions from the listed 
final states. Assuming that only one of (5 l [ R is nonzero at a time, we can then obtain constraints 
on b l [ R after specifying £ associated with a given Z 1 mass. Thus for m z , = 150 GeV, in which case 
0.008 < |f | < 0.038, 

Kr\ < °- 17 > \ h %R\ < °- 41 ' Kn\ ^ °' 44 - ( 28 ) 



z-*w 




2m z — mf 



(mj 



ml) 2 



3m| 
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For the higher masses, m z , = 0.5 - 2 TeV, we find that the upper bounds are again approximately 
proportional to their masses, 

\jf(i I ker I \y*r I 

i-^- < 1.4 , < 3.5 , < 3.8 (29) 

m z , m z , m z , 

in units of 10~ 3 GeV -1 . Stricter constraints can come from some of the other processes we study 
in the following. 



B. ll — >• 3e, r — > 3e, and r — > 3/z 



The decay /i — > e e + e receives contributions from diagrams involving the Z and Z' . From 
C int in Eq. (9), we derive the amplitude for the tree- level contributions to be 



M 



/^->-3e 



m% 



+ 



mi 



L " ;: V "/',/<'-„/> • 



gggg 
mi 



+ 



er^eeVV 



+ 



ml 



m\, J 

+ -^)eYP R ee' % P Lf , 



(e e') . 



(30) 



Here we use e and e' to distinguish the two electrons in the final state. The minus sign in the above 
equation comes from Fermi statistics. Using j3 e ^ R = t^b e ^ R and ignoring the electron mass, we can 
write the resulting branching ratio as 



B(n->3e) = 



1536 7T 3 



up 



ee 



m\ 



rrir 



UP 



ee 
^R 



°R 



rn 



+ 



mr 



+ 



mr 



+ 2 



z 
ml 



m r 



K\ 



+ 



h ee 



m. 



\ U R 



(31) 



where Tfl is the li lifetime and p e L e R = g L)R /c^ + t ? bf R . 

To evaluate the upper limits on \b e £ R \ 2 from the data on /j, — > 3e, one can try to look for 
nonzero minima of the coefficients of \b e ^ R \ 2 in the B(/jl — > 3e) formula. After scanning the 
values of £ and b e £ R satisfying the experimental requirements discussed in the previous section, 
we find for m z , = 150 GeV that the minimum of the coefficient of \b e ^\ 2 is 4.7 x 10~ 4 at 
(&l»&j?>0 - ±(0.0042,-0.021,0.025), whereas that of \b e ^\ 2 is 3.1 x 10~ 4 at {b e L e ,b%,£) ~ 
±(0.0045,-0.018,0.031). From the measured bound B{ll~ ->■ e'e+e') < 1.0 x 10~ 12 [18], we 



then extract in the m ; 



150 GeV case 



\bf I < 4.6 x 10" 



\b%\ < 5.7 x 10~ 5 . 



(32) 



For m z , = 0.5 - 2 TeV, taking similar steps we obtain the limits to be roughly proportional to m z , 
according to 



\bT 



^ < 1.4x 10~ 7 GeV" 1 , 



\bT 



Z± < 1.8 x 10~ 7 GeV" 1 . 



m z , 



m z , 



(33) 
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In the analogous case of r~ — > e~e + e~, the expression for the branching ratio can be simply 
derived from that for B(fi — > 3e) by replacing each \x in the indices with r. The same can be said 
about the coefficients of [fef^l 2 in the B{r — > 3e) formula. It follows that the measured bound 
B(t~ ->■ e-e+ e -) cxp < 2.7 x 10~ 8 [18] yields for m z , = 150 GeV 

|&L | < 0.018 , \b%\ < 0.022 , (34) 

whereas for m z , = 0.5-2 TeV 

I b eT I I b eT I 

1 L 1 < 5.3 x 10" 5 GeV" 1 , < 6.9 x 10~ 5 GeV" 1 . (35) 



As for r — > \i fi + fi , upon scanning the allowed values of b^ R and £ we find that the coefficients 

of 

^lr\ 2 m the £>(r — > 3fi) formula have minima which are vanishingly small. Consequently, this 
mode cannot provide useful restraints on \b^ R \ separately. 



C. t — y fiee and r — >■ efifi 

Another transition that can happen in our Z' scenario is r" 
bution to its amplitude is 



— > ji e + e . 
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t— >fj,ee 



b?V£ 



mi 
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rrir. 
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m|, 



eYP L efi%P L r + 



bfb R T 



It leads to the branching ratio 



rrir 



rrir 



rrir. 
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°L °R 
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The tree-level contri- 

eYP L em v P R T 



(36) 
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(37) 



where final lepton masses have been neglected and terms containing | /?c e /^c^ | = ^fl^c^el f° r 
C, C = L,R have been dropped because t| <C m z /m z ,. To determine the upper bounds on \b^ R \ 2 , 
one can again then try to seek nonvanishing minima of their coefficients in Eq. (37) which are the 
same, under the assumption that (3£ e ' eT are absent. Thus for m z , = 150 GeV we place the minimum 
to be 4.9 x 10~ 5 at (&| e ,^ e ,£) ~ ±(0.0043,-0.019,0.028). From the experimental information 
B(t~ ->■ At~e+e~) cxp < 1.8 x 10~ 8 [18], we subsequently extract for m z , = 150 GeV 



\V£ R \ < 0.019 . 



(38) 
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Similarly, for m z , = 0.5-2 TeV we arrive at 

< 6 x 1(T 5 GeV" 1 . (39) 



m z , 



Assuming b^ R = instead, we get 



m|, 



< 1.0 x 10" 8 GeV" 2 . (40) 



The constraints in the last equation are weaker by —3 orders of magnitude than those put together 
from Eqs. (32)-(35). 

For r" — >■ e~fi + fi~, the expression for the branching ratio follows from that for B(r — > jj,ee) 
with e and fx being interchanged in the indices. In this case the coefficients of |^|^| 2 in B(t — > efifx) 
have vanishingly small minima. Hence useful upper-bounds on these couplings are not available 
from B(t~ — > e~ fi + fi~) cxp < 2.7 x 10~ 8 [18]. On the other hand, assuming = we can extract 

< 1.3 x 10~ 8 GeV- 2 , (41) 
which are also very weak compared to what can be deduced from Eqs. (33) and (39). 



D. t — ^ ee/x and r — > efxfx 

Like the preceding ones, the r~ — > e~e~fi + decay receives tree-level contributions proceeding 
from Eq. (9), but involves two flavor- changing vertices exclusively. The amplitude is given by 

- (e O e') . (42) 
Neglecting the terms involving (3^(3^ as before, we consequently have 



B(r^eefi) = V |^| 2 + ^' +/' ~ ' K\ ) ■ (43) 

1536 7r d \ m z , 1 1 m|, 11 1 



The measurement i3(r — > fi + e e ) < 1.5 x 10 8 [18] then implies 



exp 

ke^t ^er I I ^ eT I 

1 L > R L > Rl < 6.8 x 10~ 9 GeV" 2 , 1 L ' R RM < 9.6 x 10~ 9 GeV" 2 . (44) 
m z , m z , 

For r~ — > e + fi~fi~, following analogous steps we obtain from £>(r~ — > e + /i~/i~) cxp < 1.7xl0~ 8 [18] 
that 

1 L ' R LM < 7.2 x 10~ 9 GeV" 2 , 1 L ' R RM < 1.0 x 10" 8 GeV" 2 . (45) 
m z , m z . 

All these results are again less strict than the corresponding constraints inferred from Eqs. (33), 
(35), and (39) by roughly 3 orders of magnitude. 
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E. Muonium-antimuonium conversion /x+e — >• /x e+ 

The experimental information on yU + e~ — >■ /i~e + is available in terms of the effective parame- 
ter G c which is defined by [18, 21] 

C cS = y/%G c Jii / P a eJi lv P c ,e + H.c. , (46) 

with C = L or R, and has been measured to be \G C \ < 0.0030 G F [18], where G F is the Fermi 
coupling constant. Attributing this to the Z' implies that 



K,r\ _ 

far less restrictive than Eq. (33). 



2\/V2\G c \ < 4.4 x 10~ 4 GeV _1 , (47) 



F. Flavor violating e+e — >■ 11' 

At e + e~ colliders, new physics could trigger the production of flavor-violating events with e/x, 
er, and /xt in the final states. In our Z' scenario, the tree-level amplitude of e + e~ — > for 
I'^l is 



m v — s 



erf (fit P L + fi R l P R )lhM e P L + ^PrV 



m z — t 

+ (Z^Z 1 , fi^b) , (48) 

where s = (j> e+ + p e -) 2 is assumed not to be close to m 2 z z , and t = (p e+ — p l+ ) 2 - The first 
experimental limits on the cross sections a (IV) = a(ee — > IV) + a(ee — > IV) were acquired by the 
OPAL Collaboration [22] at LEP-II energies, (9(200 GeV). More recent bounds on the cross sections 
at much lower energies, around 11 and 1 GeV, were reported by the BaBar [23] and SND [24] 
Collaborations, respectively. Since the theoretical cross sections tend to grow significantly as the 
energy increases from 1 to 200 GeV, the OPAL data [22] a(e/i) cxp < 22 fb, ff(er) cxp < 78 fb, and 
<7(/rr) cxp < 64 fb for the average cross sections over 200 GeV < < 209 GeV impose potentially 
stronger restraints than the others. The cross sections at these energies being more sensitive to the 
effects of m z , = 150 GeV than to those of m z , > 0.5 TeV, we discuss only the case of the former, 
in which for I — /x or r 

d(el) ~ [38.4 ((5f f € ) 2 + 175 (fi e R e t^) 2 - 20.9 (3 e L e b e L e t ? + 230 fi e R e f c 
+ 9.68 (b e L e ) 2 + 103 (b%) 2 ] \b l [\ 2 x 10 4 fb 
+ (L^R) , (49) 

v(jjlt) ~ { [193 (Pft^) 2 + QQ9/3 e L e b e L e t^ + 581 (b e L e ) 2 + (L ^ R)] |6f | 2 

- (413/3£ e t c + 717b e L e )Re(b e /;b1 T b T L e ) + (233 \bf \ 2 + 448 \b e ^\ 2 ) \b T L e \ 2 } x 10 3 fb 
+ (L R) . (50) 
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Minimizing the coefficients of |&£#| 2 ^ n °"( e an d comparing the latter to their data then yields 

|6f | < 0.76 , \b R e \ < 0.52 , \b T L e \ < 1.4 , \b T R e \ < 1.0 . (51) 

In an analogous way, (J (fir) in the absence of b e £ R gives \b^ R \ < 1.2, whereas assuming b^ R = 
instead leads to 

Wl,r\ < 0-017, |^6- L | < 0.012. (52) 
These are all weaker than their counterparts from Eqs. (32), (34), and (38) by 50 times or more. 

V. CONSTRAINTS FROM LOOP-GENERATED PROCESSES 

A. (j, — > e~f, t — > e~f, and r — ¥ fj.'j 

The flavor-violating radiative decay / — > I'-y occurs at the loop level, and its amplitude takes 
the general gauge-invariant form 

AWr = ie;kJ'(zi l P L + Z R l P R )an , (53) 

where k is the momentum of the outgoing photon, the parameters T}'^ R depend on the loop contents, 
and o UUJ = |;[7 I/ ,7 W ]- This leads to the branching ratio 

B{l^ll) = 167rm s +\Vr\)- (54) 

where r z is the I lifetime. 

This decay receives Z- and Z'- induced contributions via the diagram displayed in Fig. 3, with 
internal lepton j. Since the masses m u ,- of the external and internal leptons are small relative 
to m z z ,, it is a good approximation to retain only the lowest order terms in expanding the loop 
functions in terms of m lv -/m z z ,. In that limit, we can employ the results of Ref . [25] to derive for 
negatively charged leptons 

PlPI™>i) + (Z^Z',P^b) , (55) 




/' 



Z,Z' 



FIG. 3: Diagram for Z and Z' contributions to flavor- violating radiative decay / — > l fr y. 
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247T 2 m| 



L m V 



" /; 247r 2 m| 



where the sum is over j = e,/i,T and f3 J ^ R = (3 l i* R . Since m T ~ 17 ^> m e , we consider only 
the most enhanced terms in T, l [ l R , the ones proportional to m T . Accordingly 



e p m T b^V? 
87r 2 m|, 



e p m T bj T b T £ 
87r 2 m\, 



(56) 



Ef 



24tt 2 



e p m r 
24vr 2 



ml 



+ 



m 



TO, 



(57) 



and E^' MT follow from Yf^ tlT 



with L and R interchanged, where we have also neglected terms with 



The newest information from recent searches for these modes is B{ji —¥ e7) cxp < 2.4 x 10 
reported by the MEG Collaboration [15]. With the aid of Eqs. (54) and (56), it translates into 
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1 L > RR > Ll < 2.6 x 10" 11 



TOT 



(58) 



These numbers are 2 orders of magnitude smaller than the corresponding ones combined from 
Eqs. (35) and (39) and therefore complement them. 

The present bounds for the other 2 decays are not as strong, B(r — > e7) cxp < 3.3 x 10~ 8 and 
Bir — > /i7) cxp < 4.4 x 1CT 8 from BaBar [13, 18]. Setting b e ^ R = first, one can try to evaluate 
from these data the biggest \b e Y R T | by seeking the minima of their coefficients in the B(t — > ej, fx'j) 
formulas. Thus for m z , = 150 GeV the strongest limit we can come up with is |&| T | < 0.2, whereas 
for m z , = 0.5-2 TeV we get 



\K\ 



171 r 



< 3.3 x 10" 4 , 

< 5.9 x 10" 4 , 



\h er \ 
\°R | 

m z , 

\°R | 
to ?, 



< 5.5 x 10" 4 , 

< 7.4 x 10" 4 , 



(59) 
(60) 



all of which are less strong than the results in Eqs. (34), (35), and (39) by about an order of 
magnitude. Assuming b e l R = and b^ R = instead leads to, respectively, 



\b efl b^ T I 

1 L ' R RM < 3.6 x 10~ 7 , 



TOt 



\b^ e b eT I 

1 L ' R 2 RM < 4.2 x 10~ 7 , 



Tflr 



(61) 



which are very weak compared to the corresponding constraints deduced from Eqs. (33), (35), 
and (39) 

We remark that these / — > I'j decays effected by the Z and Z', plus additional loop-induced 
transitions I — > whose amplitudes vanish for a real photon, also contribute to the flavor- 
changing decays I — > I' I" I". However, due to the loop suppression they are less important than 
the tree-level contributions already discussed in Section IV. 
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B. Anomalous magnetic moments 



The effective Lagrangian representing the anomalous magnetic moment a l and electric dipole 
moment d t of a negatively-charged lepton I is 



Ar 7 = l 



e p a { 
Am, 



~2 ^5 ) ° l K 



(62) 



where F vu = d v A^ — d^A v is the photon field-strength tensor. Nonstandard effects of the Z and Z' 
on a l and d t appear at one-loop level, arising from the same diagram as in Fig. 3, but with /' = /. 
From Eqs. (53) and (55), we then arrive at the amplitude 



■Hi7 487r 2 m 2 



+ (Z ->• Z', £ -> 6) , 
where fc is outgoing. In view of Eq. (62), the terms without 7 5 yield 



TO; 



a, = 



12tt 2 to! 



E[ 3Re (^) m ^"(^ 



1/3 



TO; 



(63) 

+ (Z -)• Z', £ ->■ 6) . (64) 



The same expression can also be derived from Ref. [26]. Since the experimental information on a T 
is still limited [18], we will address only the / = e and fx cases. We then have from Eq. (64) 

m e m T Re(f)f !)£) ,, m„m T Re(f>f f^") 



47T 2 TO? 



47T 2 TO?, 



(65) 



where we have kept only the terms proportional to m T and also neglected terms containing /3 l £/3^. 

The SM prediction for a e agrees with its measurement, their difference being a c ^ v — af M = 
(—206 ± 770) x 10~ 14 [27]. On the other hand, the SM and experimental values of presently 
differ by about 3 sigmas, a" xp - aJ M = (29 ± 9) x 10~ 10 [27]. Consequently, we may impose 

-9.7 x 10" 12 < af < 5.6 x 10~ 12 , < ajf < 3.8 x 10" 9 , (66) 
which translate into 

-4.2 x 10- < Re( "' X) 



< 2.4 x 10~ 7 , 



Re(b IXT b TIJ ') 

< „ V L R { < 8.0 x 10~ 7 . (67) 



to 2 , GeV"' ~ ~ to 2 , GeV~ 

The result for Re(6^ r 6^ M ) is comparable to that found in Ref. [12]. These bounds are less stringent 
than those inferred from Eqs. (35) and (39), respectively. 



C. Electric dipole moments 



By comparing the 7 5 terms in Eqs. (62) and (63), the Z and Z' contributions to the electric 
dipole moment (EDM) are given by 



df = E 



e P m i 

8tt 2 



Im(/^g) Im(6g&g) 

1 * o 



TOT 



TOT 



(68) 
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Obviously, the couplings with j = I, which are real, do not matter in this case. 

Since leptonic EDM's have not yet been detected, we will again deal with only the I = e and \x 
cases, the experimental limits on d T being the least restrictive. We then have from Eq. (68) 

, _ e p m T Im(&f&g) z , _ e p m T Im(^ r 6^) 

de ~ 8n*m 2 z , ' " ~ 87r 2 m|, ' 1 ' 

where we have neglected terms containing m or f3 l i/3 J ^. Since the SM predictions 

dl M < lxl0~ 38 ecm, dJ M < 3.3xl0~ 25 ecm (70) 

are still negligible compared to the data [18] 

Id J < 1.6 x 10 _27 ecm , \d„\ < 1.8 x 1CT 19 ecm , (71) 

I e I exp — 7 I i 1 1 cxp — 7 v ' 

we can assume that the latter are saturated by Z' effects. This translates into 

KL R>l < 3.6 x 10~ 12 GeV~ 2 , l - K \ R 1 1 < 4.1 x 10~ 4 GeV~ 2 . (72) 



The first one of these was also evaluated in Ref. [12], and their result is roughly similar to ours. 
This \m{b e £b T ^) constraint appears much stricter than the one inferred from Eq. (35). But the 
comparison is actually less clear here due to the presence of a phase difference between and b T ^ 
in the former. In contrast, the \m.{b^b T ^ limit is weaker at least by 3 orders of magnitude than 
that implied by Eq. (39). 



VI. PREDICTIONS 



We summarize here the strongest limits on the Z' couplings which we have determined. Defining 
b*L% = ^L,B,l m z'i we nave f° r m z' = 150 GeV 

-4.7 x 10~ 4 < bf < 0.4 x 10~ 4 , -6.6 x 10" 4 < 1% < -0.6 x 10~ 4 , 
-2.2 x lO" 4 < bf < 5.4 x 10~ 4 , -2.0 x 10~ 4 < Vg < 6.3 x 10" 4 , 

-4.6 x 10~ 4 < b T L T < 1.6 x 10~ 4 , < b% < 5.6 x 10" 4 . (73) 

\b?\ < 3.1 x KT 7 , \b e g\ < 3.8 x 10~ 7 , 
\bf | < 1.2 x 10~ 4 , \b%\ < 1.5 x 10" 4 , 

\b^ R \ < 1.2 xlO- 4 , (74) 



while for m z , = 0.5-2 TeV 

-5.1 x 10~ 4 < b e l < -1.2 x KT 4 , -5.4 x 10~ 4 < b% < -1.1 x KT 4 

-4.3 x 10" 4 < bf < 3.4 x 10" 4 , -4.3 x 10~ 4 < Vg < 2.1 x 10" 4 , 

-6.1 x 10~ 4 < b T L T < -2.0 x 10~ 4 , 1.9 x 10~ 4 < b 7 ^ < 5.9 x 10" 4 





< 








< 







\V%\ < 1.8 x 10- 7 , 
5.3 x 10" 5 , \b%\ < 6.9 x 10~ 5 , 



(75) 



fcl $ 6 x 10- 5 , (76) 
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where all the numbers are in units of GeV -1 . The Z-pole and LEP-II measurements together have 
supplied the constraints on the flavor- conserving couplings. The numbers for the flavor-changing 
couplings have come from /i — > 3e, r — > 3e, and r — > jiee data. In addition, from // — >■ cy 

\h: R b T R " L \ < 2.6 x 10- 11 GeV- 2 , (77) 

complementary to the individual limits on b e ^ R T . Based on the results above, we now make pre- 
dictions for the largest values of a number of observables, including some of those discussed in the 
preceding two sections. Our results below can serve the purpose of guiding experimentalists in 
future searches for Z' signals. 

With these couplings, one can obviously get the decay rates of the Z' into a pair of charged 
leptons, although not their branching ratios, as we have left its couplings to other fermions unspec- 
ified. Since ^z'^vi — (|^l'| 2 + |^| 2 ) m z'/(24vr), for the flavor- conserving modes we seek values of 
the couplings which maximize the rates, but simultaneously satisfy the Z-pole and LEP-II require- 
ments discussed in Section III. For most of the Z' masses considered, the results can roughly be 
represented by 

r z ,^ e+e - < 7 x 1(T 9 m|, GeV" 2 , 
r z ,^ +M _ < 4 x 1(T 9 m|, GeV" 2 , 

T Z ,^ T+T - < 9 x 1(T 9 m\, GeV" 2 , (78) 

the exceptions being T z <^ /I + fl - < 9 x 10~ 9 m|, GeV~ 2 and T Z ,^ T + T - < 6 x 10 _9 m|, GeV -2 in the 
m Z i = 150 GeV case. For each of the flavor- violating modes, we simply choose the largest of the 
relevant set of b l ^ l R numbers in Eqs. (74) and (76) to arrive at 

IV^ e±MT < 4 x 1(T 15 m|, GeV" 2 , 
r z ,^ e±rT < 6 x 10~ 10 m|, GeV" 2 , 

r^ M±rT < 4 x 10 -10 m|, GeV" 2 (79) 

for m Z t = 150 GeV and their m z , = 0.5-2 TeV counterparts with ^ z'^T>i/ m z' ratios which are 
about 5 times smaller. 

Next are the flavor- changing Z-boson decays Z — > I'l. Since T z _^i, l ~ (|/?l'| 2 + \f3 l R \ 2 )m z /(24ii) 
and (3 l L l R ~ £ b l L l R , we again take for each mode the largest one of b l L l R from Eqs. (74) and (76), 
but employ the maximal values of £ consistent with the procedure to determine the couplings 
in Section IV. Thus, we find that the b l R l numbers for m z , = 150 GeV yield the largest branching- 
ratios, namely 

B(Z ->• e ± n T ) < 4.5 x 10 -12 , 
B(Z -> e ± r T ) < 6.8 x 10 -7 , 

B(Z -> fj^T*) < 5.1 x 10 -7 . (80) 

The latter two predictions are, respectively, only less than 25 times away from the existing limits 
B(Z ->• e ± r^) cxp < 9.8 x lO -6 and B(Z ->■ ^r^)^ < 1.2 x lO -5 [18]. 

Turning to the decays of the leptons into 3 lighter leptons, we will address only the modes that we 
did not use to derive the strictest constraints. For r — > 3/i, if the upper bounds on \b^ R \ 2 are used 
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and their coefficients in the Bir — > 3/i) expression are maximized, the resulting prediction for the 
branching ratio turn out to exceed its experimental limit. A similar situation arises in r — > ejifi, as 
can be deduced from its branching-ratio formula. Consequently, we cannot make useful predictions 
in these cases. Nevertheless, this also means that they may be potential means for probing the Z' 
within specific models. In contrast, for r — > eep, and r — > ejijj, we obtain 

B{r eep) < 1 x KT 12 , B(r ->■ e/i/i) < 7 x 10" 13 , (81) 

which come from the b^ results for m z , = 150 GeV and are much smaller than the current bounds. 
The predictions would only double if all the couplings were allowed to contribute at the same time. 
Hence the Z' effects on these 2 modes are unlikely to be detectable in the near future. 

The largest impact of the Z' on the effective coupling parametrizing the muonium-antimuonium 
conversion is also from the b£ bound for m z , = 150 GeV, 

kMe I 2 

\G C \ = < 2 x 10- 9 G F , (82) 

4v2 m z , 

far below its experimental counterpart. Accordingly, we expect that this transition is not sensitive 
to the Z' signal. 

Since the flavor-violating annihilation e + e~ — > IV depends on the center-of-mass energy, we 
will only give predictions for a(W) at 200 GeV < y/s < 209 GeV in the m z , = 150 GeV case 
to illustrate how sensitive these observables might be to the Z' signals. Thus, searching for the 
maximal rates, we get 

a(efi) < 6xl0~ 7 fb, a(er) < 0.1 fb , a(fir) < 0.05 fb . (83) 

These numbers are less than the corresponding measured bounds by about 3 orders of magnitude 
or more. 

For the radiative decays, we deal with the rates of r — )> and r — > ^7, as /i — > e*y was 
employed to produce one of the strictest constraints. Incorporating Eq. (57) in (54) and dropping 
the b e ^ R terms, we try to acquire the biggest rates by maximizing the coefficients of |^ T r T | 2 in the 
branching ratios, in a way consistent with the procedure in Section IV to extract their upper-limits, 
and subsequently applying the upper limits, one at a time. This yields 

£(r^e 7 ) < 2.3 x 10~ 8 , B{r ^7) < 2.1 x 10~ 8 , (84) 

which are close to the current limits B{r -> e7) cxp < 3.3 xlO" 8 and B{r ->■ /i7) cxp < 4.4 x 10~ 8 [18]. 

The extent of the Z' contributions to the anomalous magnetic moments and electric dipole 
moments of the electron and muon can be learned from Eqs. (67) and (72). Evidently the largest 
couplings from Eq. (74) are far from saturating the maxima of the ranges in Eq. (67) and the second 
one in Eq. (72), all drawn from comparing the SM expectations and experimental data. Since the 
first inequality in Eq. (72) involves an unknown phase difference between the couplings, nothing 
definite can be said of the Z' impact on the electron EDM in our approach. 
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VII. CONCLUSIONS 



We have considered a Z 1 boson with family-nonuniversal couplings to charged leptons and mixing 
of kinetic and mass types with the Z boson. Employing current experimental data and taking 
a model-independent apporach, we performed a comprehensive study of constraints on both flavor- 
conserving and -violating leptonic Z' couplings. Such an analysis was done for a Z' mass of 
150 GeV, as inspired by recent Tevatron anomalies, as well as higher masses of 0.5-2 TeV. We 
found that the Z-pole and LEP-II measurements together formed the strongest constraints on the 
flavor-conserving couplings. The most stringent bounds on the flavor- changing couplings came from 
the measured upper-limits of the branching ratios of the // — > 3e, r — > 3e, and r — > fiee processes. 
The radiative decay // — > supplied complementary information on the flavor- changing fi-r and 
e-r couplings. We also commented on which conditions would be relaxed if there was no mixing 
between the Z and Z'. Detailed results are summarized in the beginning of Section VI. 

With the most restricted of the extracted couplings, we computed the maximum rates of both 
flavor-conserving and -changing decays of the Z' into a pair of charged leptons as functions of 
the Z' mass. We further predicted the rates of flavor-changing Z — > W, which are not far below 
the existing measured bounds. We found that B(r — > 3/i) or B(t — > efifi) are potentially good 
observables to probe the Z' within specific models. In contrast, the rates for r — > eep, and r — > e\i\i 
were calculated to be too small to detect in the near future. Our predictions for B(r — > e^y) and 
B{r — > fij) are both very close to their current experimental limits. Finally, we commented that 
the Z' boson have relatively less impact on the anomalous magnetic moments and electric dipole 
moments of the electron and muon because of the stringent constraints on its couplings. 

Our results could also serve to constrain the rates of other Z'-mediated processes involving both 
quarks and leptons, such as the B — > X s l + l~ and B s — > decays, that have been of great 
interest recently. This would require extending the analysis to the quark sector. 

Acknowledgments 

This work was supported in part by the National Science Council of R.O.C. under Grants Nos. 
NSC-97-2112-M-008-002-MY3, NSC-100-2628-M-008-003-MY4, and NSC-99-2811-M-008-019, and 
by the National Central University Plan to Develop First-class Universities and Top-level Research 
Centers. 

Appendix A: Lagrangians with Z-Z' mixing 

The Z-Z' mixing scenario considered in this work has been described in the literature [16, 28]. 
We repeat it here using our notation for completeness. 

The interaction eigenstates for the neutral fields of the SM gauge group SU(2) L xU(l) y are, as 
usual, W 3 and B, respectively, and their coupling parameters are g and g Y . We denote the gauge 
boson of the extra Abelian group U(l)' as C and its coupling g c . Including kinetic mixing between 
B and C and mass mixing between W 3 , B, and C, we obtain the Lagrangian for the kinetic and 
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mass terms after electroweak symmetry breaking as 

1 r>voJ 



r 

*"km ~~ 



L tjuuj 75 1 /~<uu) r~i 



2 

- m^m B - m w fi W^C U + m B fi B V C V 



i« B™C UU1 + \m 2 w W£ + \m\ B 2 + \m 2 c C 2 



(Al) 



where the kinetic-mixing parameter obeys \k\ < 1 as required by the positivity of kinetic energy, 
the mass-mixing parameter jj, appears when the Higgs field carries a nonzero U(l)' charge, and 



m w = 



gv 



m B = ^ , m 2 c = M 2 +fi 2 , (A2) 

with v being the Higgs vacuum expectation value and the M c term coming from U(l)' breaking 
by a SM-singlet scalar field. Therefore, in the last equality of Eq. (Al), 



(A3) 



The kinetic part of C km can be put into diagonal and canonical form via a nonunitary transforma- 
tion: 
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(A4) 
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with 9 W being the Weinberg angle, leads to 
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one finds in terms of the mass eigenstates 

r — 1 Avw A _ l7'W7' I 1 ™2 7-2,12 7/2 

where the eigenmasses m z and m z , are already listed in Eq. (3). 

The Lagrangian for the interactions of W 3 , B, and C with fermions is 



4nt = -{9yJy 9J3 9c Jc) 



( B x \ 

V c x ) 



(A9) 



(A10) 



where J Y zc are the currents coupled to the respective fields. In terms of the fields A, Z, and Z 7 
defined in Eq. (A5), this Lagrangian can be rewritten as 



where 



— e P ^cm ~~ 9z Jz Z\ — 9z> Jz 1 Z'\ , 



int 



(All) 



e 7 A 

"em 



c w9y-Jy + s w9J^ 1 9zJz — C w9-J3 s w9y-Jy > 



jx _ 9c J c , 
9z' J z' — y 9y j y 5 



(A12) 



with 



e p ~ 9 s w ~ 9y c w > 
sinx 



9z = — » 



*x = 



cosx 



c,„ = COS 



W > 



s„, = sin 9 



w > 



sinx = k , c = cosx = Vl — k 2 



(A13) 



In Eq. (5) we reproduce only the part of C' int involving Z and Z'. We note that the field A x coupled 
to the electromagnetic current J A m is massless, as Eq. (A6) indicates, and hence identical to the 
physical photon. 



Appendix B: Cross sections of e+e — > l + l 



From the amplitude in Eq. (23), with each of the propagators now assumed to have a simple 
Breit-Wigner form, follows the cross section 
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and the forward-backward asymmetry 

-4fb 



a FB (e+e ->• l+l ) 



a(e+e- -)• l+l~) ' 

where a = e 2 /(4:7i) is the fine-structure constant, T zz , are the total widths, and 



(B2) 



a FB (e + e~ 



a 



(s — m|) 2 + r|m : 



+ 



(s - m|,) 2 + r|,m|, 



[(^f) 2 -(^) 2 ][W) 2 - W) 2 ]^ , [(&r) 2 - fe) 2 ] [(^) 2 - (b l k) 2 ]s 

^ . r / o \ o -no o "i 



64tt[(.< 



2 + r|m|] 



647r[(s-m|,) 2 + r|,m|,] 



+ 



(gggf - gfg) - Plb ll R ) (s - m|) (g - m z ,)s 

32n[(s -m 2 z ) 2 + T|m|] [(s - m 2 z ,f + r|,m|,] 



(B3) 



the lepton masses having been neglected. These formulas agree with those in the literature [29]. 
Here f3 e £^ = g L R) ^c^+t^b £ £^. In our numerical computation away from s ~ rn 2 z z ,, we set V z z , = 0. 
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